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Abstract 

Stellar fundamental properties—masses, radii, effective temperatures—can be extracted from observations of eclipsing binary 
systems with remarkable precision, often better than 2%. Such precise measurements afford us the opportunity to confront the 
validity of basic predictions of stellar evolution theory, such as the mass-radius relationship. A brief historical overview of 
confrontations between stellar models and data from eclipsing binaries is given, highlighting key results and physical insight 
that have led directly to our present understanding. The current paradigm that standard stellar evolution theory is insufficient 
to describe the most basic relation, that of a star’s mass to its radius, along the main sequence is then described. Departures 
of theoretical expectations from empirical data, however, provide a rich opportunity to explore various physical solutions, 
improving our understanding of important stellar astrophysical processes. 


1 Introduction 

Eclipsing binary (EB) systems are some of the best tools 
for testing stellar evolution theory. Advancements in stellar 
evolution theory motivated by confrontations between stel¬ 
lar model predictions and observations of EBs are second, 
perhaps, only to studies of globular cluster color-magnitude 
diagrams. The power of EB observations to further stellar 
evolution theory is derived from observer’s ability to extract 
accurate stellar masses and radii of the individual compo¬ 
nents. Being that mass is the primary input parameter in 
stellar models and radius is a primary output quantity, EBs 
permit direct comparison with model predictions at the most 
fundamental level. Furthermore, the presence of two, pre¬ 
sumably coeval, stars in a binary system forces stellar mod¬ 
els to predict the properties of both stars at the same age, 
effectively making each binary system a miniature star clus¬ 
ter. 

This review will attempt to characterize advances in low- 
mass stellar evolution theory that were and are motivated 
by investigations of EBs. Naturally, the definition of a low- 
mass star is strongly dependent on context, so for reference, 
I will refer to low-mass stars as stars with masses below 
0.8 Mg. Stars in the range of 0.8 to 1.2 M 0 will be re¬ 
ferred to as “solar-like” stars. These definitions are largely 
set by the characteristics of stellar evolution models in the 
given mass regimes. The lower boundary of 0.8 Mg is de¬ 
fined by the growing importance of non-ideal contributions 
to the gas equation of state and a need for detailed radia¬ 
tive transfer models to prescribe the surface boundary con¬ 
dition. Additionally, models below the 0.8 M 0 boundary are 
less sensitive to model input parameters such as metallicity, 
the convective mixing length, and age compared to models 
above that threshold. The upper boundary to the “solar-like” 
regime at 1.2 M 0 corresponds, roughly, to the mass thresh¬ 
old above which main-sequence stars are believed to develop 
a convective core. 

Definitions outlined in the previous paragraph are the re¬ 
sult of knowledge accumulated through decades of research 


on stellar evolution. The importance of many of the above 
effects were revealed over time, some through the study of 
EBs. It is therefore instructive to begin by placing our mod¬ 
ern understanding of low-mass stellar evolution theory into 
perspective with a recapitulation of historical progress. Fol¬ 
lowing the historical evolution of stellar evolution theory, a 
summary is provided of where we are today, what is “state 
of the art,” and where EBs are making significant contribu¬ 
tions to the advancement of stellar evolution models. Finally, 
ideas about how EBs may contribute significantly to stellar 
evolution theory are presented. 

Before embarking on the historical perspective, one must 
caution that the historical review presented below attempts 
to cover the most important advances in the state of low- 
mass stellar evolution theory. As such, individual stud¬ 
ies contributing to the the knowledge and inspiration of 
those studies cited below may appear under appreciated, as 
might the effort that goes into measuring stellar masses and 
radii. There is no intent to minimize the meticulous work 
of theorists and observers who put tireless hours into ad¬ 
vancing physical models and extracting exquisite measure¬ 
ments from their data. Without these efforts, theorists would 
be unable to construct advanced computational models and 
they would be unable to casually remark on the agreement 
between observations and theory. 

2 Historical perspective 

The history of testing low-mass stellar evolution theory 
with EBs begins with the first observations of YY Gemino- 
rum (hereafter YY Gem, also Castor C; Adams & Joy 1920) 
and its subsequent classification as an EB (van Gent 1926; 
Joy & Sanford 1926). YY Gem consists of two MO-type stars 
that are nearly identical when it comes to their fundamental 
properties ( M ss 0.60 M 0 , R « 0.62 R 0 , T e ff« 3800 K). 
Although it would be another 20 years before YY Gem pro¬ 
vided clues about the inadequacy of stellar model physics, it 
will be continually mentioned throughout the historical de¬ 
velopment of low-mass stellar evolution theory. 
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The first major contribution provided by YY Gem (and 
EBs) to low-mass stellar evolution theory came in the early 
1950s. At this point, significant advancements in energy 
generation processes via nuclear reactions had been made 
and the proton-proton (p-p) chain was established as the 
dominant energy production mechanism in late-type stars 
(Aller 1950; Salpeter 1952). However, comparisons of stel¬ 
lar evolution models computed assuming energy produc¬ 
tion via p-p chain predicted stellar luminosities far in ex¬ 
cess of observed luminosities for stars below the Sun in the 
HR-diagram (Aller 1950; Aller et al. 1952). Stromgren (1952) 
identified numerous possible shortcomings in stellar model 
physics, including the validity of the equation of state (EOS) 
and adopted radiative opacities. Most importantly, Strom¬ 
gren noted that the prevailing assumption of radiative equi¬ 
librium (e.g., Eddington 1926) throughout stellar interiors 
was likely incorrect and should be abandoned. Instead, it 
was argued that the properties of a surface hydrogen con¬ 
vection zone in late-type stars should be explored. 

Following Stromgren’s suggestion, Osterbrock (1953) 
computed the first set of stellar models of late-type stars with 
an extended outer convective envelope. Osterbrock found 
satisfactory agreement between model calculations with a 
deep convective envelope and observations of YY Gem for 
models where the hydrogen convection zone comprised the 
outer 30% of the star, by radius. The study by Osterbrock, 
motivated by observational properties of YY Gem, revealed 
that stellar models must account for the transport of energy 
by convection, a point that seems rather trivial from our 
modern point of view, but represented a significant advance 
in stellar evolution theory, at the time. However, while the 
inclusion of an outer convection zone provided agreement 
between model predictions and the properties of YY Gem, 
there was still significant disagreement between models and 
the properties of stars of later spectral type. Most notable 
were disagreements with the stars in Kruger 60, a visual bi¬ 
nary consisting to two mid-to-late M-dwarfs (today classi¬ 
fied as M2.0 and M4.0). Although the system is not eclipsing, 
its role in stellar structure theory is important and worth 
mentioning. Both stars were observed have a luminosity 
well below the predictions of stellar evolution models, much 
the same as the stars in YY Gem before the introduction of 
a hydrogen convection zone. These remaining errors were 
initially attributed to missing effects of electron conduction 
as an energy transport mechanism (Osterbrock 1953). 

Carefully ruling out severe observational errors in the 
determination of stellar properties for Kruger 60, Limber 
(1958a) pointed out that several additional pieces of physics 
were likely necessary to reconcile model predictions with 
the observations. These physics included partial electron 
degeneracy of the stellar interior and a more rigorous treat¬ 
ment of radiative and convective temperature gradients for 
computation of energy transport in stellar interiors. How¬ 
ever, even when all of these effects were included in models, 
significant disagreements with Kruger 60 remained (Limber 
1958b). Instead, Limber extrapolated on the advances pro¬ 
vided by Osterbrock (1953) and allowed the model interiors 
to be in full convective equilibrium. Under this assumption, 
models were able to provide predictions consistent with the 
properties of the stars in Kruger 60 (Limber 1958a,b). 

Precise measurement of the fundamental properties of the 
stars in YY Gem helped initiate the establishment of the cur¬ 
rent paradigm that outer convection zones grow deeper to¬ 


ward later spectral types, effectively overturning the pre¬ 
vailing hypothesis of radiative equilibrium. The revolution 
was so complete that stellar interiors of the latest type stars 
were determined to be fully convective! However, agree¬ 
ment between models and data was foreseen to be short¬ 
lived. Observational data still possessed large error bars that 
greatly helped to mask the identification of further disagree¬ 
ment with model predictions (Limber 1958b). At the same 
time. Limber warned that measuring stellar radii would lead 
to more robust comparisons with stellar evolution models. 
To this end, he suggested further advances in theory that 
may be required, including a larger nuclear reaction network 
to include the burning of lithium and deuterium, departures 
from hydrostatic equilibrium, “violent atmospheric activity” 
now associated with magnetic activity, and the influences of 
rotation and magnetic fields. Yet, despite the obvious com¬ 
plications that could enter into the picture, Limber urged 
caution (Limber 1958a, pg. 368), 

We should not attempt to introduce these added 
complications unless and until the simpler models 
can no longer account for the observations. 

Building on the framework laid out Limber (1958b) and 
following the suggestion that improvements to simpler mod¬ 
els be explored thoroughly, the next significant advance¬ 
ments in low-mass stellar evolution theory came with the 
introduction of the mixing length theory (MLT) of convec¬ 
tion for stellar models (Bohm-Vitense 1958; Henyey et al. 
1965) and improved numerical schemes for solving the set 
of stellar structure equations (Henyey et al. 1964). These ad¬ 
vances permitted more detailed stellar models to be con¬ 
structed. Several groups created sets of low-mass stellar 
evolution models that include the aforementioned advances, 
as well as more detailed stellar atmosphere calculations fol¬ 
lowing a gray T(r) construction (e.g., Copeland et al. 1970; 
Hoxie 1970). These calculations included treatment of var¬ 
ious atomic ionization states and molecular dissociation on 
the EOS and used improved opacity data for bound-free and 
free-free absorption. 

Stellar model calculations were found to be in good agree¬ 
ment with the observational mass-luminosity relationship, 
but model radii were found to under-predict the radii of stars 
by up to 40% below about 0.8 Mg (Hoxie 1970). Solace was 
found by noting that observational errors were still domi¬ 
nant, masking any disagreements (Hoxie 1973). However, 
models also failed to fit the observationally determined prop¬ 
erties of YY Gem within 3 a of the observed radius uncer¬ 
tainties, a fact that was not explicitly mentioned. Instead, YY 
Gem was used to highlight observational uncertainties in the 
single star effective temperature scale, justifying the large 
uncertainties in the empirical mass-radius relationship. Nev¬ 
ertheless, potential deficiencies in the mass-radius, and to a 
lesser extent the mass-luminosity planes motivated Hoxie to 
identify the largest sources of uncertainty in low-mass stel¬ 
lar models. Two specific pieces of physics were the inclusion 
of non-ideal effects in the EOS, which would likely affect the 
stellar radius predictions, and development of more sophis¬ 
ticated boundary conditions, such as non-gray model atmo¬ 
spheres, which would have the largest effect on the stellar 
T e g and luminosity (Hoxie 1970, 1973). 

The next twenty years saw a significant increase in the 
number of well-characterized EBs (e.g., Popper 1984). Un¬ 
fortunately, only one low-mass EB was added to the list with 
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Figure 1: Mass-radius diagram for known eclipsing binaries with masses and radii determined to better than 3% on the MS (filled points; 
Torres et al. 2010; Kraus et al. 2011; Feiden & Chaboyer 2012a) and 5% on the pre-MS (open points; Stassun et al. 2014). The theoretical ZAMS 
from Dartmouth stellar evolution models is shown as a dashed line. 


YY Gem. That particular system would, however, become 
just as crucial of a system as YY Gem for comparisons with 
stellar evolution theory. CM Draconis was characterized by 
Lacy (1977) and was found to contain two very-low-mass 
stars with masses similar to those in Kruger 60, meaning 
they were very likely fully convective throughout their in¬ 
terior (Limber 1958b), although this was not immediately 
recognized. The stars in CM Dra also appeared to he well 
above the theoretical zero-age-main-sequence (ZAMS) for 
Population I objects in a mass-radius diagram, but this was 
attributed to the fact that CM Dra may be a Population II 
object (Lacy 1977). 

The surge of well-characterized EBs during these inter¬ 
vening years were compiled by Andersen (1991). In that re¬ 
view, systems were required to have component masses and 
radii measured with better than 2% precision so that they 
might act as strong tests of stellar evolution theory. How¬ 
ever, there was only one low-mass system in that collec¬ 
tion: YY Gem (Leung & Schneider 1978; Andersen 1991). CM 
Dra’s components had masses determined to 4% and radii to 
3% (Lacy 1977) and therefore did not merit inclusion in An¬ 
dersen’s compilation. 

Low-mass stellar models made their next big leap with 
the Lyon stellar models (Baraffe et al. 1995). What set their 
models apart from other groups developing models concur¬ 
rently (e.g., Dorman et al. 1989; Burrows et al. 1993) was in¬ 
clusion of an advanced EOS designed specifically to accu¬ 
rately model the cool, dense plasma characteristic of low- 
mass stars (Saumon et al. 1995) and adoption of non-gray 
model atmospheres (Allard & Hauschildt 1995) to define sur¬ 
face boundary conditions for their interior models. The lat¬ 
ter feature, in particular, allowed their models to accurately 
treat energy transfer in optically thin regions of the stellar 
atmosphere, where convection and radiation both contribute 
significantly to the overall energy flux (Dorman et al. 1989; 
Burrows et al. 1993; Allard & Hauschildt 1995; Baraffe et al. 


1995). Naturally, one of the first tests of their models was 
comparing model predictions to the properties of YY Gem 
and CM Dra (Chabrier & Baraffe 1995). Initial results were 
encouraging, showing that their models were able to repro¬ 
duce the observational properties of both systems. However, 
revisions to the fundamental properties of YY Gem and CM 
Dra would again show models unable to accurately predict 
stellar properties (Metcalfe et al. 1996; Torres & Ribas 2002). 
Nevertheless, the physics advances implemented in the Lyon 
models (Chabrier & Baraffe 1997; Baraffe et al. 1998) remain 
state of the art, making them still highly relevant. 

3 State of the art 

In the years since Andersen (1991), there was another ex¬ 
plosion in the number of well-characterized EBs, especially 
in the low-mass regime. To synthesize the wealth of obser¬ 
vational data, an updated compilation of EBs with precisely 
known properties was published by Torres et al. (2010). That 
review contains 95 systems, but only five systems possess 
at least one component mass below 0.8 M 0 , and only one 
system contains a star below the fully convective boundary. 
The lack of low-mass systems was relieved in subsequent 
years thanks largely to photometric surveys searching for 
exoplanets and variable stars, such as Kepler (Carter et al. 
2011; Doyle et al. 2011), MOTESS-GNAT (Kraus et al. 2011), 
MEarth (Irwin et al. 2009, 2011), All Sky Automated Survey 
(ASAS; Helminiak & Konacki 2011; Helminiak et al. 2011, 
2012, 2014), Super Wide Angle Search for Planets (Super- 
WASP; Triaud et al. 2013; Gomez Maqueo Chew et al. 2014), 
and HATNet (Zhou et al. 2014). Figure 1 shows a subset of 
known EBs (see John Southworth’s DEBCat 1 for a complete 
listing) whose masses and radii are determined with preci¬ 
sions better than 3% and 5% for main sequence (MS) and 
pre-MS EBs, respectively. 


1 http://www.astro.keele.ac.uk/ jkt/debcat/ 
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Figure 2: Stellar evolution models fail to reproduce the proper¬ 
ties of low-mass main-sequence stars in eclipsing binaries (adapted 
from Feiden & Chaboyer 2012a). 

Standard stellar evolution models of low-mass stars em¬ 
ploy similar physics to the Lyon series: surface boundary 
conditions defined by non-gray model atmospheres, modern 
radiative opacity computations, and an EOS that accounts 
for partial electron degeneracy, Coulomb interactions, and 
other non-ideal effects. The rapid influx of high quality data 
and high precision mass and radius measurements permits 
a more reliable evaluation of stellar evolution model accu¬ 
racy. Such evaluations indicate that the problems identified 
by Hoxie (1973) still remain, though at a lower level than 
originally suggested. Canonically, model radii are quoted to 
underpredict observed values by between 5% and 10% (Ribas 
2006). These values can be decreased to below 5%, if po¬ 
tential metallicity and age variations are taken into account 
(Feiden & Chaboyer 2012a), as shown in Figure 2. However, 
it is worth noting that the most significant outliers are also 
the most well-characterized systems, including YY Gem (at 
0.6 Mq) and CM Dra (near 0.23 M 0 ). 

At the same time, it is clear that problems facing stellar 
evolution theory below 0.8 M 0 are not unique to this mass 
regime. Williams (this volume), for example, illustrates nu¬ 
merous problems that exist with high mass models. Per¬ 
haps more worrying, is that studies of EBs with solar-like 
components reveal that models in the solar-mass regime fail 
to reproduce observations (Popper 1997; Clausen et al. 2009; 
Vos et al. 2012). Stellar evolution theorists have previously 
comforted themselves knowing that, while there are prob¬ 
lems in various locations of the HR Diagram, models of the 
Sun and solar-like stars are accurate. If stars in EBs are repre¬ 
sentative of the single star population, these errors can have 
a profound impact on interpretation of asteroseismic data, 
on characterizations of exoplanet host stars, and potentially 
on isochrone fitting to globular and open cluster data. Given 
that modern stellar evolution models show systematic dis¬ 
agreements with EB data, what possible physics may be in¬ 
complete or entirely missing from stellar models? 

3.1 Observational evidence 

The first clue was presented by Mullan & MacDonald 
(2001), who noted that low-mass EBs showed similar depar¬ 
tures from standard stellar model isochrones as single ac¬ 
tive M-dwarfs, a conclusion supported by additional com¬ 
parisons of EBs to the single star population (Morales et al. 
2008). This lead Mullan & MacDonald (2001) to the hypoth¬ 
esize that magnetic activity may be the culprit for the ob¬ 
served model-observation disagreements. Their hypothe¬ 
sis was supported by the fact that most EBs had short or¬ 
bital periods, typically less than three days (due to observa¬ 
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tional biases; Ribas 2006). Nonetheless, stars in these EBs 
are expected to have rotational periods synchronized with 
the orbital period, meaning these stars are rapidly rotating 
(Zahn 1977). Rapid rotation then drives strong magnetic 
fields through a hydrodynamic dynamo. 

Further evidence in support of the magnetic hypothesis 
was provided when it was shown that radius discrepancies 
between models and observations of low-mass EBs correlate 
with magnetic activity in the form of coronal X-ray emis¬ 
sion (Lopez-Morales 2007). More active stars, identified as 
having a higher ratio of X-ray luminosity to bolometric lu¬ 
minosity, were found to show larger disagreements with 1 
Gyr, solar metallicity Lyon models. In that same study, there 
appeared to be no correlation between stellar metallicity and 
radius deviations for EB systems, a trend that was apparent 
for single stars and was proposed as an alternative hypoth¬ 
esis to magnetic activity (Berger et al. 2006; Lopez-Morales 
2007). This is still some of the strongest evidence in support 
of the magnetic hypothesis. 

A second approach to assessing the influence of magnetic 
activity is to directly compare model-observation radius dis¬ 
agreements with stellar rotation. Since activity is theorized 
to increase with increasing stellar rotational velocity, one 
expects to see shorter period binaries display larger radius 
deviations, while longer period binaries should show bet¬ 
ter agreement with model predictions. Kraus et al. (2011) 
found that short period binaries, with orbital periods less 
than 1 day, displayed a significantly higher level of radius in¬ 
flation compared to model predictions than those with peri¬ 
ods greater than 1 day. This was supported, in part, by a later 
study that found a significant break in the level of radius in¬ 
flation at 1.5 days (Feiden & Chaboyer 2012a). However, no 
significant difference was found for EBs with orbital peri¬ 
ods of 1 or 2 days, suggesting the difference may be due to 
low number statistics, particularly for systems with orbital 
periods greater than 3 days. Even accounting for the poten¬ 
tial influence of rotation on convection, through the use of 
the Rossby number, 2 no significant correlation between ra¬ 
dius deviations and rotation is observed (Feiden & Chaboyer 
2012a). Curiously, Feiden & Chaboyer (2012a) still uncover a 
correlation with levels of coronal X-ray emission, similar to 
Lopez-Morales (2007). However, their correlation is neither 
as clear nor as strong as in Lopez-Morales (2007), and it is 
observed to be driven largely by two data points, including 
YY Gem. 

Presently, there appears to be no clear observational con¬ 
sensus identifying either magnetic activity or rotation as 
drivers of the observed radius discrepancies. A few binary 
systems with orbital periods greater than 18 days are now 
known whose stars show significant radius disagreements 
with stellar models (see, e.g., Feiden & Chaboyer 2013a). 
This does not immediately rule out magnetic activity and/or 
rotation as a culprit, but certainly adds a layer of complexity 
to problem. These stars in long period EBs, which show ra¬ 
dius disagreement, tend to have masses placing them below 
the fully convective boundary, perhaps suggesting rotation 
and magnetic fields are driving errors for models of stars 
with radiative cores and convective outer envelopes and that 
models of fully convective stars are plagued by other errors. 
This evidence has, in any case, inspired theoretical investi- 


2 Ro = Prot/iconv, the ratio of the convective turnover time to the ro¬ 
tational period 
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gations into the effects of magnetic fields on low-mass stellar 
structure. 

3.2 Theoretical advancements 

The theoretical basis behind suspecting magnetic fields is 
that strong magnetic fields can suppress the flow of energy 
across a given surface within a star (e.g., Thompson 1951; 
Gough & Tayler 1966). Since stars must maintain a given 
flux at their surface, commensurate with the energy produc¬ 
tion rate in their core, the inhibition of the outward trans¬ 
fer of flux at any point within the star will force the star to 
readjust to compensate for the suppressed flux. This com¬ 
pensation is expected to take the form of an inflated radius 
and decreased T e ff. The same physical argument applies to 
both strong magnetic fields globally suppressing convection 
(e.g., Mullan & MacDonald 2001) or the blocking of radiative 
flux near the surface by starspots covering some fraction of 
the stellar surface (Spruit 1982a). 

3.2.1 Magnetic suppression of convection 

Stellar evolution models that include parametrized de¬ 
scriptions of the influence of global magnetic fields have 
been constructed by Lydon & Sofia (1995), Ventura et al. 
(1998), D’Antona et al. (2000), Mullan & MacDonald (2001), 
and more recently by Feiden & Chaboyer (2012b). In¬ 
vestigation of the effects on low-mass stellar properties 
in the context of model disagreements with fundamen¬ 
tal stellar properties were carried out with the mod¬ 
els by Mullan & MacDonald (2001) and Feiden & Chaboyer 
(2012b), the latter of which uses the prescription de¬ 
rived by Lydon & Sofia (1995). Results from both groups 
suggest that magnetic suppression of convection is able 
to reconcile model radius predictions with observations 
of low-mass benchmark EBs, most notably YY Gem 
(Feiden & Chaboyer 2013a; MacDonald & Mullan 2014) and 
CM Dra (MacDonald & Mullan 2012; Feiden & Chaboyer 
2014). 

In stars with radiative cores and convective envelopes, 
such as YY Gem, quantitative predictions of surface mag¬ 
netic field strengths necessary to reconcile model radii ap¬ 
pear consistent with empirical data. Model magnetic field 
strengths range from a few hundred gauss up to a few kilo- 
gauss. Direct measurements of magnetic field strengths 
on active, early M-dwarfs confirm these predictions are 
realistic (Reiners 2012), as do estimates of magnetic field 
strengths from indirect indictors, such as chromospheric Ca 
ii emission and coronal X-ray emission (Feiden & Chaboyer 
2013a; MacDonald & Mullan 2014). Interior magnetic field 
strengths predicted by models (~ 10 - 100 kG) also ap¬ 
pear consistent with expectations from more realistic mag¬ 
netohydrodynamic simulations of stellar magnetic fields 
(Brown et al. 2010). However, models of stars with radia¬ 
tive cores are relatively insensitive to deep interior magnetic 
field strengths, as the influence of magnetic fields in the su- 
peradiabatic layer near the stellar surface is of greater con¬ 
sequence (DAntona et al. 2000). The broad consistency of 
model magnetic field strengths with empirical data suggests 
that magnetic models may be capturing relevant physics for 
these types of stars, despite simplified magnetic field pre¬ 
scriptions. 

The situation for models of fully convective stars, typ¬ 
ified by CM Dra, is more complex. Surface mag¬ 


netic field strengths predicted by models span a simi¬ 
lar order of magnitude, anywhere between about 0.5 kG 
to 5.0 kG (MacDonald & Mullan 2012; Feiden & Chaboyer 
2014). These values are reasonable when compared to aver¬ 
age surface magnetic field strengths of mid-to-late field M- 
dwarfs (~ 3 kG; e.g., Reiners & Basri 2007). Debate about 
whether magnetic fields are actively inflating fully convec¬ 
tive stars largely focuses on predicted interior magnetic 
field strengths (Chabrier et al. 2007; MacDonald & Mullan 
2012; Feiden & Chaboyer 2014). Unlike stars with radia¬ 
tive cores, models of fully convective stars are less sensi¬ 
tive to the strength of the magnetic field in the near-surface 
layers and appear to be relatively more influenced by in¬ 
terior magnetic field strengths (Mullan & MacDonald 2001; 
MacDonald & Mullan 2012; Feiden & Chaboyer 2014). Inte¬ 
rior magnetic field strengths typically need to be in excess 
of 1.0 MG. Interested readers are encouraged to consult the 
aforementioned references for a full account of the argu¬ 
ments for and against the appearance and maintenance of 
megagauss magnetic fields in fully convective stars. It is 
likely that stronger constraints from magnetohydrodynami- 
cal simulations are needed to further the debate. Neverthe¬ 
less, there is some concern over the magnitude of interior 
magnetic field strengths and, thus, global magnetic fields as 
a solution for the mass-radius disagreements between mod¬ 
els and observations. 

3.2.2 Magnetic activity - starspots 

Alternatively, it may not be global magnetic fields inhibit¬ 
ing convective flows that suppresses flux and forces stars to 
inflate, but intense local concentrations of magnetic fields 
producing starspots on the stellar photosphere. Magneti¬ 
cally active stars show modulations in their lightcurves due 
to the presence of dark spots. 3 Dark spots reduce the radia¬ 
tive output from the stellar surface by trapping excess en¬ 
ergy at their base (Spruit 1982b). Depending on the lifetime 
of a given spot and the efficiency at which the trapped en¬ 
ergy can be redistributed, a star may grow larger and cooler 
in response to the presence of spots (Spruit 1982a). 

It was shown that, in the case of low-mass EBs, effects 
due to spots may be more significant in driving structural 
changes in a stellar model than effects from global inhibition 
of convection by magnetic fields (Chabrier et al. 2007). This 
is of particular consequence for models of fully convective 
stars, like CM Dra, where global inhibition of convection 
may not be sufficient, at least with realistic magnetic field 
strengths (Chabrier et al. 2007; Morales et al. 2010). How¬ 
ever, starspots may have an additional influence on stud¬ 
ies of EBs. The presence of starspots can bias radius mea¬ 
surements from EB lightcurves (largely through biasing of 
the measured radius sum) toward larger radii by up to 6%, 
depending on the properties of spots (Morales et al. 2010; 
Windmiller et al. 2010). For systems like YY Gem and CM 
Dra, Morales et al. (2010) estimated that the stellar radii may 
be overestimated by about 3%. It is not possible to necessar¬ 
ily attribute all of the observed errors to this bias, but assum¬ 
ing these stars have smaller radii assists theoretical models 
of radius inflation through magneto-convection or starspot 
flux suppression by lowering the required amount of radius 


3 There is the distinct possibility that lightcurve modulations are caused 
by bright spots, but it has been argued that this is unlikely (Berdyugina 
2005). 
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inflation (Morales et al. 2010; MacDonald & Mullan 2012). 

Critical to the idea that spots bias radius measurements is 
that spots are preferentially located near the rotational poles 
and cover a significant fraction of the stellar surface (be¬ 
tween 40 - 60%). Suppression of flux by spots at the surface 
is not sensitive to distribution, but is highly dependent on 
surface coverage and spot temperature contrasts. Presently, 
there is no strong empirical evidence for polar spots on ei¬ 
ther rapidly rotating low-mass stars or fully convective stars 
(see Section 4.3 in Feiden & Chaboyer 2014, for an in-depth 
discussion). In fact, there is evidence that spots must be 
randomly distributed across the surface of fully convective 
stars to produce the diversity of lightcurve morphologies for 
stars in the young open cluster NGC 2516 (Jackson & Jeffries 
2013). This leads to the possibility that fully convective stars 
have random distributions, but large filling factors (total sur¬ 
face coverage). Observations of FeH molecular band features 
in spectra M-dwarfs reveal that M-dwarfs are probably cov¬ 
ered nearly everywhere by 1 kG magnetic fields, with small 
patches of intense 5 kG - 8 kG magnetic fields, averaging 
to about 3 kG over the surface (Shulyak et al. 2011), sup¬ 
porting this idea. However, this assumes a one-to-one cor¬ 
relation between the presence of 1 kG magnetic fields and 
the appearance of spots. Investigations of solar magnetic 
fields reveal that 1 kG magnetic fields may not be sufficient 
to inhibit convective flows (Mathew et al. 2004). Convective 
flows in low-mass stars likely react differently to the pres¬ 
ence of magnetic fields than convective flows in the Sun, but 
it raises the question of whether starspot properties required 
by stellar models are realistic. 

4 Exploring alternative solutions 

While magnetic fields and activity are able to provide a 
solution to the mass-radius problem with low-mass stars, 
questions remain about the reality of stellar evolution model 
predictions. These questions will undoubtedly be answered 
by future observations. However, perhaps an equally valid 
epistemological approach is to explore alternative theoreti¬ 
cal solutions with the aim of removing the need to invoke 
magnetic fields, or with the aim of ruling out the other alter¬ 
natives and thus bolstering arguments in favor of the mag¬ 
netic hypothesis. Recalling discussions in Section 2, this is 
the approach advocated by Limber (1958a). 

4.1 Convection 

One alternative to the magnetic hypothesis is that con¬ 
vection in low-mass stars is significantly less efficient than 
it is in more solar-like stars. This argument has been ex¬ 
plored previously, but investigations have largely focused on 
individual systems. Suggestions for why convective prop¬ 
erties may be different among individual EBs that show 
radius disagreements with models include magnetic fields 
(Cox et al. 1981), rotation (Coriolis force acting on convec¬ 
tive flows; Chabrier et al. 2007), and intrinsic differences 
due explicitly to stellar properties (i.e., mass, metallicity, 
T e g; Lastennet et al. 2003). While it is important to demon¬ 
strate that manipulating convective properties—here largely 
the convective MLT parameter cumlt - provides relief to 
noted model-observation disagreements, it makes identi¬ 
fying physical explanations for those manipulations diffi¬ 
cult. Instead, statistical properties from a sample of EBs pro¬ 
vides an opportunity to reveal meaningful trends that betray 


physics associated with changes in stellar convective prop¬ 
erties. 

Such a study was recently performed for solar-like EBs 
(Fernandes et al. 2012), where cimlt was manipulated to 
provide the best agreement between models and EBs. Re¬ 
sults show that q;mlt does not appear to be mass depen¬ 
dent, as stars of equal mass from different EB systems require 
different values for «mlt- Instead, Fernandes et al. (2012) 
demonstrate that omlt correlates with v sin i, where faster 
rotating stars require lower ckmlt values. While the results 
are tantalizing, comparison of models against solar-like EBs 
introduces errors due to unknown stellar ages and compo¬ 
sitions (both helium abundance Y and metal abundance Z), 
which must be simultaneously fit. Although stars in a given 
EB must lie along the same isochrone, which isochrone is 
assumed correct can decidedly alter modeling conclusions. 
Extending this type of study to low-mass EBs would, how¬ 
ever, help relieve modeling uncertainties associated with 
unknown stellar ages and compositions, as low-mass stel¬ 
lar models are relatively less sensitive to assumptions about 
these parameters. Having at least some constraint on stellar 
metallicity can drastically improve the quality of model com¬ 
parisons, so observers are anyway encouraged to put forth 
effort to measure stellar metal content. 

It is quite remarkable that, through a systematic compari¬ 
son of the low-mass EB population to stellar models, it may 
be possible to extract information regarding the hydrody¬ 
namic properties of low-mass stars as a function of a range 
of variables, such as mass, [M/H], T e g, and usin* in a sim¬ 
ilar manner as is done with asteroseismology (Bonaca et al. 
2012). This knowledge would lead to significant improve¬ 
ments in the treatment of convection in stellar evolution 
models and would represent yet another facet of stellar 
physics advanced by studies of EBs. 

4.2 Heavy element composition 

Stellar metallicity plays only a minor role in modern dis¬ 
cussions of the low-mass EB mass-radius problem. How¬ 
ever, it is not clear, at least to the author, that metallicity ef¬ 
fects are not partially responsible for the observed disagree¬ 
ments. Although Lopez-Morales (2007) failed to identify a 
correlation between metallicity and radius inflation among 
low-mass stars in binaries, the lack of a trend was largely 
driven by the points associated with CM Dra. Revisions to 
the fundamental properties of CM Dra (Morales et al. 2009) 
and recent estimates of its metallicity (Terrien et al. 2012; 
Kuznetsov et al. 2012) have led to a shift in the location of 
CM Dra in the metallicty-radius inflation diagram, revealing 
circumstantial evidence that stellar metallicity may be re¬ 
lated to observed radius errors among fully convective stars 
(Feiden & Chaboyer 2013b). 

However, correlations with metallicity among fully con¬ 
vective stars in EBs does not follow the same pattern as was 
observed for single stars (Berger et al. 2006). Instead of more 
metal rich stars showing larger radius inflation, radius errors 
among fully convective stars in EBs appear to negatively cor¬ 
related with stellar metallicity, with more metal poor stars 
displaying greater radius errors. This may point toward er¬ 
rors in the model EOS or opacities. However, the sample 
is small, with only six stars in four systems, and is largely 
driven by the location of CM Dra in the diagram. The no¬ 
tion is bolstered slightly by the fact that it does not matter 
what metallicity is adopted for CM Dra, it still falls along 
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the observed negative correlation. It should be noted that a 
recent investigation by Zhou et al. (2014) claims to find no 
trend of radius inflation with metallicity among fully con¬ 
vective stars in EBs. However, their most precisely measured 
star, and their most metal-poor ([M/H] = —0.6 dex), falls 
along the relation suggested by Feiden & Chaboyer (2013b). 
Confirmation of the properties of other stars in their sample 
is required. As the number of fully convective stars in EBs 
with metallicity estimates is increased, it will become clear 
whether hints of the trend are spurious or not. 

4.3 Helium abundance 

An interesting possibility is to use low-mass EBs to in¬ 
vestigate stellar helium abundances (Limber 1958b). Helium 
abundances in low-mass stars cannot be measured directly, 
and therefore specification of Y in stellar models relies on 
prior assumptions regarding the relationship between Y and 
Z. The significant role that helium plays in governing the 
structure and evolution of low-mass stars, coupled with our 
ignorance of Y in low-mass stars, makes it a noteworthy 
suspect in our attempts to solve the mass-radius problem 
(Valcarce et al. 2013). Although observational confirmation 
of model predicted Y values cannot be obtained, consistency 
with other measures of local helium abundances in the solar 
neighborhood can be sought. 

As with studies of adjusting o:mlt> previous investiga¬ 
tions largely focused on individual systems, such as CM 
Dra (Paczynski & Sienkiewicz 1984; Metcalfe et al. 1996) and 
UV Psc (Lastennet et al. 2003). Populations of solar Tike EBs 
were used to probe the helium enrichment as a function 
of metallicity in the solar neighborhood (Ribas et al. 2000; 
Fernandes et al. 2012). It is usually assumed that helium 
abundance is linearly proportional to metallicity, such thats 

Y(Z) = H' Z, (1) 

where Yp is the primordial helium abundance. Solar-like 
EBs suggest that the slope of the relation is AY/AZ = 
2 ± 1 (Ribas et al. 2000; Fernandes et al. 2012), consistent 
with estimates from single K-dwarfs in the solar neighbor¬ 
hood (e.g., Casagrande et al. 2007). However, the inferred 
primordial helium abundance from these studies is well be¬ 
low the primordial abundance estimated from Big Bang Nu¬ 
cleosynthesis (Yp ss 0.225 compared to Yp, bbn ~ 0.249; 
Peimbert et al. 2007). 

Inferring helium abundances involves simultaneously fit¬ 
ting multiple stellar model parameters (i.e., cumlt, Y, Z, age) 
to provide the best possible agreement between model pre¬ 
dictions and observed properties of stars in EBs. It was men¬ 
tioned in Section 4.1 that these parameters have a strong im¬ 
pact on stellar model calculations of solar-like stars, intro¬ 
ducing degeneracies in the optimization problem. Presently, 
very few EBs have metallicity estimates, making it difficult to 
add observational priors in the optimization scheme. Metal¬ 
licity estimates would help tremendously, particularly for 
solar-like stars. At the same time, with the recent increase 
in the number of well-characterized low-mass EBs, it is now 
possible to use low-mass EBs to constrain Y(Z). Models 
are considerably less sensitive to fit parameters, particularly 
age, as compared to solar-like stars. Metallicities are still 
a concern, but methods of estimating bulk metallicities for 
low-mass stars are showing promise (Rojas-Ayala et al. 2012; 
Mann et al. 2013). Assuming how helium abundances scales 


with metallicity is a crucial in stellar evolution modeling, 
and studies of low-mass EBs provide the best chance to re¬ 
veal this relation. 

4.4 Probing internal structure 

Perhaps the ultimate test of stellar evolution theory that 
EBs can provide is a direct inference of a star’s internal den¬ 
sity structure through the measurement of apsidal motion in 
EBs with eccentric orbits. This is very fitting for this occa¬ 
sion, where we gather to remember the contributions to bi¬ 
nary star science by Zdenek Kopal. Observing the precision 
of an EB orbit’s periastron with sufficient accuracy and pre¬ 
cision yields an average interior structure constant, k 2 of the 
two component stars. For any individual star, &2 quantifies 
the object’s central mass concentration (Kopal 1978). Objects 
that are very centrally condensed, with tenuous outer layers 
and dense cores compared to the average density (i.e., high 
mass stars) are described by lower values than objects that 
have a more even distribution of mass (i.e., M-dwarfs). 

The fact that observations only reveal the component av¬ 
eraged /c 2 value means that equal mass binaries are prefer¬ 
able. Fortunately, low-mass binaries are more likely to be of 
equal mass (Bate, this volume). However, low-mass stars 
have a lower fraction of binarity and typically form with 
smaller semi-major axes (see Bate, this volume), meaning 
they will tend to rapidly circularize from tidal interactions 
(Zahn 1977). All considered, the properties of low-mass bi¬ 
naries makes the chances of finding suitable systems for de¬ 
riving &2 quite low. 

Despite this, if apsidal motion in a low-mass binary is ac¬ 
curately measured, we will learn a great deal about the ac¬ 
curacy of stellar models. Primarily, we will learn whether 
models predict accurate density profiles for the deep interi¬ 
ors of low-mass stars. The apsidal motion constant is largely 
determined by the deep interior structure, with smaller in¬ 
fluences from the near surface layers, where the density in 
individual layers is much less than the average density of the 
star. If models accurately predict k '2 for stars where models 
do not predict accurate fundamental properties (radius, T e fj), 
then the problems can be isolated to the near surface layers. 
In contrast, if the deeper layers in models are found to be 
inadequately described, this points towards a different set 
of physics, such as opacities and the EOS. At least one low- 
mass binary gives us this opportunity, KOI-126 (Carter et al. 
2011), where the apsidal motion constants may be deter¬ 
mined for two fully convective stars with a precision of 1%. 

5 Conclusion 

Our current understanding of modeling errors are very 
much like the early view on M-dwarf stars (Russell 1917, 
commenting on Kruger 60), 

It is obvious that these stars exhibit every charac¬ 
teristic which might be supposed typical of bodies 
at the very end of their evolutionary history, and on 
the verge of extinction. 

At the time, this view of M-dwarfs was well justified, but 
nonetheless it now appears rather strange, as M-dwarfs are 
thought to live for hundreds of billions of years, meaning 
they are actually at the very beginning of their evolution¬ 
ary history. Indeed, it is very easy to grow confident that 
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observed errors between models and stars in EBs are them¬ 
selves becoming extinct. However, as history of comparing 
models to theory has continually demonstrated, more often 
than not, errors persist even with the most sophisticated 
models. It is entirely likely that, as M-dwarfs are actually 
at the very beginning of their lives, we are only at the be¬ 
ginning of our journey of reconciling errors between stellar 
evolution theory and observations. It is clear, though, that 
well-characterized EBs, like YY Gem, will continue to play 
critical roles in advancing stellar evolution theory. 
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